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I. INTRODUCTION 

In the last two decades since the seminal papers by Berezinski and by Kosterlitz and 
Thouless, the critical behavior of the two-dimensional XY (or plane rotator) model has 
been much studied numerically, mainly on the square lattice [Q-|^ 0^13' ^^^' sometimes, 
also on the triangular lattice 0J^. The steady increase of the computers power, important 
recent progress in MonteCarlo (MC) algorithms |l^ and the calculation of 0(/5^°) high tem- 



perature expansions (HTEs) on the square lattice [§], have produced increasingly accurate 
verifications of the Kosterlitz and Thouless (KT) theory [J15|,|l6|. However, since on rea- 
sonable grounds it has been disputed @,0 that these studies are really conclusive, further 



quantitative evidence is still valuable. 

Here we present an extension (by two terms up to order /?^^) and a new analysis, by the 
methods of Ref. [§, of HTEs for the XY model on the triangular lattice. Our study gives 
further support to the KT picture and leads to fairly precise estimates of the KT critical 
parameters. On the contrary our results appear not to be consistent with a conventional 
power-law critical behavior unless the exponents are very large, namely 7 ^ 3.2 and 7-|-2z/ ^ 
5.4. If this is the case, our 0(/5^^) series do not seem to be sufficiently long to yeld accurate 
estimates of the critical exponents. 

In the second Section, we recall briefly the main predictions of the KT theory. The third 
Section is devoted to an analysis of the series by ratio-extrapolation, and by rational and 
differential approximants techniques. The last Section contains our conclusions. 

II. THE MODEL AND THE HT SERIES 

The Hamiltonian of the two-dimensional XY model is 

H{s} = -^-J:s{x)-s{x'). (1) 

{x,x') 

Here s{x) is a two-component classical spin of unit length associated to the site with position 
vector a; of a 2-dimensional triangular lattice and the sum extends over all nearest neighbor 



lattice pairs {x,x'). 

Our series have been computed by a FORTRAN code which solves recursively the 
Schwinger-Dyson equations for the correlation functions [[7|,p!7|. 



Here we shall analyze the HTE of the spherical moments of the correlation function 
m^'-^P) for / = ( reduced susceptibility ) and / = 2. 

Table I shows the HTE coefficients through /?^^ of the nearest-neighbor spin-spin cor- 
relation function. In Tables II and HI we have reported the expansion coefficients for the 
moments m^^\l3) with / = and 2. 

The main predictions of the non-rigorous renormalization group analysis of the plane 
rotator model [p!5| , p!6[ can be summarized as follows. 

The correlation length ^(/3) is expected to diverge as /3 | /?c with the unusual singularity 

m^U{P)=exp{^)[l + 0{T)] (2) 

where t = [3c — (3. 

The value of the exponent a predicted in Ref. ITH] is a = 1/2 and 6 is a non- universal 
positive constant. 

At the critical temperature, the asymptotic behavior of the two-spin correlation function 

as r = |x| ^ oo is expected to be 

< s(0) ■ s{x) >oc -^^^^^^[1 + 0(ln(ln(r))/ln(r))] (3) 

The values predicted [|1^,0 for rj and are, respectively, rj = 1/4 , 9 = 1/16. 

From Eqs. (0) and @, together with the usual scaling ansatz, it follows that, for / > r]—2, 
the correlation moment m'^'-' (/?) should diverge as /? t /3c with the singularity 

m«(/5) oc r~'U{Pf-'^'[l + 0(rhn(r))] (4) 

At Pc a line of critical points should begin which extends to /5 = oo, so that for /3 > /3c 
both ^ and the correlation moments remain infinite. 

Finally it is worth mentioning the rigorous bound [^ 



/3c > 2/?f = hn3 = 0.5493... (5) 

where /3^ is the inverse critical temperature of the Ising model on the triangular lattice. 



III. AN ANALYSIS OF THE HT SERIES 



We estimate the critical parameters by simple modifications of current methods of series 



analysis 0j|,|g^ 



In Refs. d^J^ we have shown that the ratios rn{m^^'^) = a^'V'^n+i of the successive HTE 
coefficients of the correlation moment m^'-* {(3) for large n behave as: 

r„(m«)=/3,+ ^^^ + 0(l/n^) (6) 

with e = YT^, Ci = — ((2 — r/ + l)abj3cY and A = mm(l, 2e). According to the KT prediction 
we should have e = 2/3. This is a clear signature of the KT singularity and it should be 
detected in a ratio analysis of the HTEs. 

On the other hand if, instead of (0) and (Q), we had conventional power-law critical 
behavior so that, as /5 t /3c, 

m(')(/3) ~ r-T-''^[A + Bit"^ + ..] (7) 

where, A > and, as usual, 7 and u denote the susceptibility and the correlation length 
exponents respectively, we would obtain a formula analogous to Eq. (|^) with e = 1 and 
A = 1 + A, namely 

r„(m«) =(3^+ ^-^^'^-^'^^ + Oil/n'^^) (8) 

In Fig.l we have plotted versus 1/n the sequences of ratios r„(m'^')) for / = and 2. 
These ratio plots exhibit a residual curvature and an increasing slope for large n. If Eq. (H) 
were a correct representation of the asymptotic behavior of rn{m^''^), the 0{l/n) terms in 
(I) should be suppressed by forming the linearly extrapolated sequences 

r^i)(m(')) = nrn{m^^^) - (n - l)r„_i(m(')) = [3^ + 0(1/^^+^) (9) 



which, for large n, should approach with vanishing slope their common limit j3c- This does 
not happen (at least up to n = 14), as it is shown in Fig.l where we have also plotted the 
sequences r^^)(m'-')) versus 1/n. The estimates of /5c obtained from r^\m^^^) still increase 
rapidly with order. 

Turning to critical exponents, we have computed a sequence of (unbiased) estimates of 
7 + /z/ by the formula 

_ (n - l)Vn(m(')) - n{n - 2)rn_i(m('^) 
^^^^^"~ nr„_i(m(0)-(n-l)r„(m(0) ^^^ 

The sequences of estimates so obtained for 7 and 7 + 2z/ are plotted versus 1/n in Fig.2. 
We conclude that, under the assumption of power-law scaling, the simplest extrapolations 
suggest 7 > 3.2 and 7 + 2z/ > 5.4. These estimates for the exponents are larger than those 
obtained from a fit of ( square lattice) MC data to power law critical behavior. 

Let us assume now that Eq. (|^) is valid instead of Eq. (^, then by reporting the r„(m'-')) 
sequences versus 1/ra^/^, we should obtain nicely straight plots, as indeed is observed in Fig. 3. 
We can suppress the 0{l/n^^^) terms in the sequences rn{m^^^) by forming the (nonlinearly) 
extrapolated sequences 

Unfortunately the sequences obtained are not regular enough to give a much more precise 
estimate of (3c by a further extrapolation in 1/n. The results are reported in Fig. 3. We can 
infer that (3^ = 0.683 ± 0.004. 

A direct unbiased estimate of e in terms of ratios is obtained from the sequence 

en = nH^^) (12) 

where tn = Vf) • If the ratios rn{x) and r„(x^) have the asymptotic behavior (^, the 
sequence e„ will provide estimates of e. A quantity m„ analogous to tn may be defined in 
terms of the moment m*^^)(/5) and its square, and, via Eq. |T2|, the corresponding sequence 
e^ may be formed. The sequences e„ and e^ have been plotted versus 1/n in Fig.4. We 



have also reported the sequence e^ as computed from the susceptibihty xi of the triangular 
lattice Ising model [^, in order to emphasize the qualitatively different behavior of the two 
cases. These tests support the KT theory and are inconsistent with a power-law singularity: 
if that were the case, as it appears from the Ising model plot, the limit of the sequences 
should be 1. An extrapolation of e„ and ej^ to n = oo by the Barber-Hamer method pO 



leads to the estimate a = 0.51 ± 0.04, while for the sequence e„ we get a = 0.02 ± 0.04. 
Notice that this test is able to distinguish sharply power-like from KT singularities, unless 
asymptotic behaviors have not yet set in. This is true even if the critical exponents are large: 
for instance a reasonable model series, like xl ( having 7 = 3.5), produces an e-sequence 
with the same qualitative behavior as e„. 

The rest of our analysis uses differential approximants (DAs) or simply Pade approxi- 
mants (PAs). 

The expected singularity structure of /n(x), namely ln(x) = ^{P)iPc — P)~'^ + B(P) with 
A{P) regular and -B(/5) at most weakly singular at /3c, should be reproduced with reasonable 
accuracy by inhomogeneous first order DAs. We have selected DAs [n/l; m] with 1 < n < 6, 
2 < m < 5 and 2 < / < 6. From the non-defective approximants of the reduced sample 
which use at least 12 series coefficients, we get the unbiased estimates Pc = 0.680 ± 0.002 
and a = 0.49 ± 0.03. If we set a = |, we get the biased estimate p^ = 0.681 ± 0.002. 
Conservatively we have estimated uncertainties as three times the standard deviation of 
the sample. The results remain essentially unchanged whether we use the Fisher- Au Yang- 



Hunter-Baker or the Guttmann- Joyce definitions pO[ of the DAs. 

We have also computed the PAs to the logarithmic derivative of ln(x)//5. This quantity 
should discriminate between the structures (^ and (|^) of the critical singularity, since the 
residues at the critical poles have either to approach a, if (^ holds, or to vanish, if (^ holds. 
From the PA table for the location of the critical pole of the approximants to Dln[ln(x)//3] 
and the PA table for the residues we get the estimates jSc = 0.684±0.003 and a = 0.54±0.04. 
If we set cr = i, we get the biased estimate Pc = 0.680±0.003. The PA and the DA estimates 
are therefore consistent, the small difference in the central values being probably due to 



background effects. 

Due to the slower convergence of the m^^-* series a similar DA analysis of the correlation 
length (using ln(4-) ) is unsuccessful unless we assume a = |. In this case we get the 
estimate (3c = 0.684 it 0.004 consistent with the previous ones from x, but somewhat less 
accurate. 

Assuming a = |, the non- universal parameter b may be obtained by computing PAs of 
the quantity 

C{(3) = ^irrHl + m(2)/x) = b + Oir^ (13) 

at f3 = pc- Taking (3^ in the range [0.680, 0.682], we estimate b = 1.27 ± 0.05. 

Finally, the critical index t] governing the large distance behavior of spin-spin correlation 
function may be estimated (without bias on a) by PAs of the ratio 

ln(x) 2-r] 

at /3 = jSc- Allowing as above for the uncertainty on j3c, we estimate rj = 0.27 ± 0.01. 

If we assume a power-law singularity (|^), from a study of the location of the singularities 
of the PAs to the logarithmic derivative of the susceptibility Dln{x), we should be able to 
estimate (3c, and from their residues, the critical exponent 7. As we have already observed 
in the case of the HTEs for the square lattice p|, both the PA tables for the poles and for 
the residues contain many "defective entries" or blanks and show a very poor convergence. 
These features of the PAs suggest that the critical singularity is not a power or, at least, 
that our series are still too short. If we insist in producing anyway some estimate of the 
critical parameters, then, by averaging over all relevant entries of the PA tables for the poles 
and residues of the approximants to D\n{x) we get (3c = 0.655 ± 0.008 and 7 = 3.7 ± 0.6. 
This estimate for 7 is consistent with those from ratio tests, but larger than those obtained 
in Ref. M from power-law fits to MC data. 



IV. CONCLUSIONS 

Let us finally compare our results to those obtained in previous studies of the model on 
the triangular lattice. To our knowledge, no MC simulation is available for the ferromagnetic 
XY model on the triangular lattice although, by virtue of the higher coordination number 
of this lattice, the approach to scaling behavior is expected to be faster than in the square 
lattice case, for the same lattice size. On the same grounds one can also argue |^ that, for 
a given number of HT coefficients, the triangular lattice series are "effectively longer" than 
the square lattice ones. In fact, we expect from our analysis that it would perhaps take only 
an 0(/?^^) triangular lattice series to reach the same precision as with our 0(/5^°) square 
lattice series ||8||. 



An extensive review of the square lattice numerical studies can be found in Refs. ||8|JTl 
and needs no duplication here. As to HT studies ||l|] on the triangular lattice, we recall that 
the early ones, based on ten term series, were essentially inconclusive and did not provide 
reasonably stable estimates of the critical parameters. When the HT series were extended 
to twelve terms [^, an analysis by the four-fit method gave convincing indications of the 
validity of the KT predictions and yelded the estimates /?c = 0.687 ± 0.009, a = 0.5 ± 0.1 
and r/ = 0.27 ±0.03. 

We believe that our extended series and our new tests further substantiate the KT 
picture and provide more precise estimates of the critical parameters, otherwise consistent 
with those of Ref. p. We have pointed out that the study of ratio plots, and of the PAs 
to the logarithmic derivative of x ^^^ ^, are strongly suggestive that the critical singularity 
is not a power. If, however, we want to pursue a power-law interpretation, we should not 
overlook the unusual facts that, for a closely packed lattice, 0(/5^^) series do not seem long 
enough to determine with reasonable accuracy the critical parameters 7 and u and that the 
rough estimates obtained for these exponents are significantly larger than those observed in 
conventional critical phenomena and show a trend to increase with the number of coefficients 
used. These estimates are also larger than those obtained in the power-law fit of Ref. to 

8 



square lattice MC data. However we have reexamined the most recent and extensive MC 



simulation data [y,|Tl|] on the square lattice and we have observed that they are consistent 
with a power-law finite size scaling Ansatz if we choose the larger exponents suggested by 
our series ( it seems that this range of parameters has not been explored in the fits of Refs. 
PJTT[|). Our interpretation of all these facts has been already pointed out at the conclusion 
of our analysis of the square lattice series p|: as the critical point is approached, larger and 
larger effective exponents are needed to fit an infinite order singularity by an ordinary power 
singularity. 

We may conclude that, if we extrapolate the HT series consistently with the KT behavior, 
we find stable values of the critical parameters in good agreement with the KT fits to the 
MC data and with our previous studies of HTEs on square lattice. Our unbiased estimates 
of the critical parameters are (3^ = 0.680 ± 0.002, a = 0.49 ± 0.03 and rj = 0.27 ± 0.01. 

If we set a = ^ we obtain the estimate Pc = 0.681 ± 0.002 and by also fixing Pc at its 
central value we find b = 1.27 ± 0.05. 
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FIGURES 
FIG. 1. The successive ratios of the HTE coefficients of various moments are plotted versus 1/n. 

The ratios r„(x) are represented by sohd squares; r„(m*-^)) by sohd triangles. We have also plotted 

the linearly extrapolated ratio sequences rn (x) (open squares), and rn [mP') (open triangles). 

FIG. 2. Unbiased estimates of the critical exponent 7 of the susceptibility under the assump- 
tion of a power-law critical singularity obtained from the ratios r„(x) (open squares). Analogous 
estimates of the exponent 7 -|- 2z^ as obtained from r„(m^^') (open triangles). 

FIG. 3. Ratio plots for the HTE coefficients versus 1/n^'^. The ratios r„(x) are represented 
by solid squares; rn{m^^') by solid triangles. The ratio sequences have been extrapolated in l/n?'^ 
obtaining the sequences s„(x) (open squares), s„(m(^)) (open triangles). A further extrapolation 
in \/n of the sequences s„ gives sii {x) (open circles), s„ (m^^^) (crosses). 

FIG. 4. The sequences e„ (open squares), e^ (open triangles), as computed from the quantities 
tn introduced in Eq. |l^ , and from the analogous ones u„, are plotted versus 1/n. The dashed line 
indicates the KT prediction for the value of e. We have also reported the analogous sequence e^ 
(crosses) as computed from the susceptibility of the triangular lattice Ising model. 
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TABLES 





TABLE L HTE coefficients of the nearest neighbor correlation function. 


order 


coefficient 



1 0.500000000000000000000000000000 

2 0.500000000000000000000000000000 

3 0.437500000000000000000000000000 

4 0.312500000000000000000000000000 

5 0. 182291666666666666666666666667 

6 0.072916666666666666666666666667 

7 -0.039550781250000000000000000000 

8 -0.204481336805555555555555555556 

9 -0.448027886284722222222222222222 

10 -0.763769531250000000000000000000 

11 -1.139343883373119212962962962962 

12 -1.581509003815827546296296296296 

13 -2.114185611785404265873015873016 

14 -2.760278320043385554453262786596 



14 



TABLE II. HTE coefficients of the susceptibility m^^' . 



order coefficient 

1.000000000000000000000000 

1 3.000000000000000000000000 

2 7.500000000000000000000000 

3 16.87500000000000000000000 

4 35.62500000000000000000000 

5 72.06250000000000000000000 

6 141.2734375000000000000000 

7 270.1728515625000000000000 

8 506.3834635416666666666666 

9 933.5703776041666666666666 

10 1697.512101236979166666666 

11 3050.264278496636284722222 

12 5424.862119886610243055555 

13 9561.162654477074032738095 

14 16716.55156094636866655299 
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TABLE III. HTE coefficients of the second correlation moment m^'^'. 

order coefficient 

0.000000000000000000000000000 

1 3.000000000000000000000000000 

2 18.00000000000000000000000000 

3 72.37500000000000000000000000 

4 239.6250000000000000000000000 

5 703.7500000000000000000000000 

6 1902.406250000000000000000000 

7 4835.969726562500000000000000 

8 11719.07975260416666666666666 

9 27326.64085286458333333333333 

10 61726.45278320312500000000000 

11 135743.4609174940321180555555 

12 291741.0660864935980902777777 

13 614640.8839340452163938492063 

14 1272465.830598210733403604497 
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